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Abstract 

We study the effects of a scalar condensate on a class of 2+1 dimensional non- 
Fermi liquids by introducing fermionic probes in the corresponding asymptotically 
AdS^ black hole backgrounds. For the range of parameters and couplings we consider 
we find gapless fermionic excitations whose properties are model-dependent. 
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1 Introduction and Conclusions 



An important problem in condensed matter physics is the understanding of fermionic 
systems at strong coupling. The strongly coupled nature of these systems, and nu- 
merical instabilities due to the fermion sign problem make it hard to investigate such 
systems using conventional analytic or numerical methods. A non-Fermi liquid picture 
is believed to emerge in some situations; in these cases the fermionic system has a well- 
defined Fermi surface but differ significantly in other respects from the Landau Fermi 
liquid. Non-Fermi liquids are interesting for a number of reasons. For example, it has 
been proposed that the normal state of high-T c superconductors, and metals close to 
a quantum critical points are examples of non- Fermi liquids [DI2]- 

In string theory, one can use the gauge-gravity duality to investigate various strongly 
coupled systems, which map a strongly coupled field theory to a weakly coupled grav- 
itational system [3]. Recently there has been much progress in finding gravitational 
analogues of condensed matter and fluid dynamics phenomena, see for example [3HH] 
and references therein. 

An interesting recent development is the work by Liu et al [9], where the authors 
have studied probe fermions in an extremal AdS black hole geometrjjf). The authors 
calculate the retarded Green function of the fermionic field, showing that the system 
behaves like a non-Fermi liquid, or in other words in a certain range of parameters a 
well-defined Fermi surface seems to exist. Some interesting non-Fermi scaling behavior 
can be demonstrated [9|[T5] . 

The main motivation of this note is to study of the effects of a scalar condensate on 
the fermionic syste To this end, we extend the black hole background by turning on 
a scalar condensate; In this paper we will mostly be concerned with the non-extremal 
geometries. Since the condensate strength can become large, it becomes necessary to 
consider its effects on the gravitational background. We next introduce probe fermions 
in this geometry with a suitable coupling to the scalar. The coupling to the scalar, and 
other parameters we consider here (such as the fermion charge and the temperature) 
extend and generalize previous work [15]. We find that the system always contains 
sharp fermionic excitations, which are gappless and become stable in the zero temper- 
ature limit. The detailed properties of those fermionic excitations, depend strongly 
on the parameters chosen, in a way we discuss in detail. The robustness of the Fermi 
surface is perhaps surprising, and might be related to the large N limit. 

The plan of this note is as follows. We describe the required background in sections 
2 and 3. Section 2 is devoted to introducing the models of holographic superconduc- 
tivity we will use here. Section 3 introduces the tools needed to probe those systems 
using bulk fermions. We then introduce features of the Fermi surface (as exhibited 
by properties of the fermion spectral functions) is sections 4 and 5: in section 4 we 
discuss the high temperature case, and in section 4 we discuss the low temperature 
(condensed) phase. 

In this note we concentrated on the physics of the holographic superconductors ob- 
tained from the simplest Abelian Higgs model in the bulk [21H24] , There is another set 
of models utilizing non- Abelian gauge fields in the bulk [25], which are more interesting 
in some ways. For instance, the bulk couplings are more constrained, and the model 

6 For other work on holographic non- Fermi liquids see [imimriTMl4| ■ 
7 For previous work on such effects see [15H20] . 
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allows for both s-wave and p-wave |26| superconductor^. We hope to report on our 
results studying this set of models shortly [29] . 

2 Asymptotically AdS^ Black Holes 

The set of backgrounds we are interested in asymptote to an AdS^ geometry, with 
possibly some profiles for the gauge and scalar fields. This corresponds to a 2+1 
dimensional dual field theory at a finite temperature, with a finite chemical potential, 
and possibly resulting in a condensation of a scalar operator. We are interested in 
various such geometries, which dominate the thermodynamics of the dual field theory 
at different temperature ranges. 

The system is described by gravity coupled to a Maxwell field and a charged scalar 
field (an Abelian Higgs model), with the Lagrangian: 

C = R + J2 ~ \ F ^ F ^ ~ I W " mM\ 2 ~ V(\^\), (2.1) 

where F is the electromagnetic field strength, and qb is the charge of the scalar field 
ip. We follow the conventions in [24] . 

Assuming solutions with spherical symmetry, the metric has the general form: 

ds 2 = -g(r)e- x( - r) dt 2 + -^r + r 2 (dx 2 + dy 2 ) (2.2) 

9{r) 

with A = (j){r)dt and ip = ip(r). By a suitable gauge choice we can assume ip to be real, 
and the scaling symmetries of the metric allow us to set L = 1. The full non-linear 
equations of motion are given by: 

*'' + (£-* + ?>' + ^-^ = o 

5'+(i-^) 5 + ^-3r + ^ = 0, (2.3) 

with the simple potential V(ip) = ^m^ifj 2 where is the mass of the scalar field. 
Furthermore, we will mostly take the scalar ip to be massless (m^ = 0), unless otherwise 
indicated. 

The above field equations have two scaling symmetries that will turn out to be 
useful. They are: 

r -)■ ar, (t,x,y) ->• (t,x,y)/a, g -> a 2 g, <j> -> a</>, 

^ -> a 2 e*, t at, <f> -»• <f>/a. (2.4) 

We are interested in non-extremal black hole solutions. We use the first symmetry to 
set the black hole horizon radius r + to 1, and the second symmetry to set \ to zero at 

For recent work on this set of models see [27l[28] . 
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(a) Scalar ip (b) Metric function g 

Figure 1: The profiles for g(r) and ip(r) at T e ff = 0.036. For reference, the grey monotoni- 
cally increasing curve is for the pure RN black hole. 



the boundary. With these choices, the following boundary conditions can be used to 
fix the remaining free parameters in Eq. (|2.3j) : 



0(r+) = 0, <f/(r+) = E, g(r+) = 0, 
V>(r+) = V+, tf(r+) = —, 1 V '^ +) i „ r (2.5) 

With the above choice of the boundary conditions, the solutions are determined by two 
parameters: ip+ and E. The latter quantity can be interpreted as the electric field at 
the horizon. From the equations of motion Eqs. (12. 3D . the general asymptotic form of 
the scalar field ip near the boundary is of the form: 

The parameters 0+ and E in Eqs. (|2.5p can be traded for ^i, ^2- In what follows we 
choose to set "01 to zerc(^|. This fixes E, and the remaining free parameter ip + can be 
used to vary the temperature of the solution. 

In Fig Q] we plot typical profile for the scalar fied ^p(r) and the function g{r). The 
temperature is T e fj = ^- = 0.036. The expected asymptotic forms for this solution 
are ip ~ 1/r 2 and g ~ r 2 for large r. 

As the temperature is reduced, the thermodynamics is dominated by one of two 
different backgrounds. Above the critical temperature, there is no scalar hair, and 
the relevant background is the Reissner-Nordstrom (RN) black hole. Below the critical 
temperature, the scalar begins to condense; though around the critical temperature this 
can be treated as a small perturbation, as we lower the temperature further backreac- 
tion ultimately becomes important. Thus we have to use a non-extremal backreacted 
black hole with scalar hair, constructed in [241 . As the temperature approaches zero 



9 We set either ipi or ip2 has to be zero, as we are interested in the field theory in the absence of external 
sources. Massless scalars allow for two inequivalent quantization, where either ipi or ip2 is interpreted as a 
normalizable mode. We choose the quantization for which -02 is normalizable; the alternative quantization 
gives results which are qualitatively similar to the ones discussed below. 
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the system resembles the zero temperature backreacted black hole with a zero radius 
horizon, constructed in [30]. We find that in the zero temperature limit our results 
approach those of [31], as they should. 



3 Fermions 

We now concentrate on the effect of a scalar condensate on non-Fermi liquid behavior. 
The basic idea, introduced in [9], is to introduce fermions in our geometries, whose 
Green's functions probe the existence and properties of a Fermi surface in the boundary 
theory. In our case there can be also a scalar field present in the bulk, and therefore 
we have to specify the coupling between the fermions and the scalar condensate. 
We take the action for the bulk fermion ^ to be: 

= J d d+1 x^g~ i(^T M V M y - mf * - A|V>| 2 **), (3.7) 

where T M are the curved space Gamma matrices and: 

* = Vm = Q M + ^ i0abM r ab - iq f A M . (3.8) 

Here M denotes bulk spacetime indices while a, b denote tangent space indices. Greek 
letters denote indices along the boundary directions. Thus r ab are the tangent space 
Gamma matrices. We also choose units such that R = 1 in the AdS geometry. The 
charge of the fermionic field is denoted by qf, and the fermion mass is m. We are 
mostly working with massless fermions, m = 0. 

The last term in Eq. (j3.7p is a quartic coupling between the scalar ip and the Dirac 
fermion and A is a tunable parameter controlling the coupling to the scalar, as- 
sumed to be positive. This is the most general coupling for general choices of scalar 
and fermion charges, though cubic couplings can exist when q b = 2qj |31| . The quar- 
tic coupling can be absorbed into an effective (radial dependent) mass term for the 
fermions: 

M(r) = m + A|V>| 2 (3.9) 

Since the quartic coupling does not necessitate a definite ratio of fermion and boson 
charge, we are free to vary their ratio. In the following we keep the charge of the 
complex scalar field to be q b = 1, which is in the range where a phase transition 
occurs, but keep the freedom to vary qj. 

With the quartic coupling only, the mechanics of the fermionic equations remains 
more or less the same |9j. We choose the following basis for the Gamma matrices and 
the (4 component) spinors: 

where ^± are two-component spinors and 7 M are (2+l)-dimensional gamma matrices. 
We can now separate the radial and boundary coordinate dependencies in ^ as follows: 

*± = (-gg rr )~h- iwt+ik > xl <S>±, $± = ft^) (3.11) 
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For the 7 M , we choose the basis 7 = io^T 1 = ci,7 2 = 03- We also use the rotational 
symmetry of the system to set &2 = 0. The Dirac equations then reduce to two sets of 
decoupled equations: 




— {d r =F M^fg~r~r~)y± = Ti{h - u)z T , 



' — {d r ± M*fg r ~ r ~)z T = ±i{ki + u)y±, 

9rr 



with 



u 



V -9tt 



(3.12) 



(3.13) 



We define the ratios = iy~/z + , 
be written as: 



'9. 



—drt± = -2M^e± t (fci t «) ± (fci ± «)e±- 

The retarded Green's function Gij is given in terms of the quantities £± by: 

-2.U A - 



iz^/y + , in terms of which eqs. (13.120 can 

(3.14) 



Gb = lim e 



^ 



Gn 
G 22 



(3.15) 



The spinors f± satisfy infalling boundary conditions at the black hole horizon, which 
is located at r + = 1. This implies: 



lim £± (r) 



r— >i 



(3.16) 



This completes our review of the setup. We will be interested in the properties of the 
Fermi surface as we vary the temperature, or the quartic coupling a 1c ] In order to do 
that, we solve numerically the system of equations (|3.14p in the various backgrounds 
obtained by solving the bosonic equations, Eqs. (|2.3p . From this, we can calculate the 
retarded Green's function (Eq. (|3,15p ) and the fermionic spectral densities (related 
to the imaginary part of the retarded Green's function), and investigate the effects of 
turning on A and varying the temperature. 



4 Reissner-Nordstrom black hole (T > T c ) 

The Reissner-Nordstrom black hole exists for any temperature, and is the dominant 
phase for T >T C . This background was originally discussed in (9j, and we concentrate 
here on the behavior at non-zero temperature, to set up our notation and provide a 
comparison to other backgrounds. The Reissner-Nordstrom metric is given by: 

ds 2 = -f{r) dt 2 + ^—+ r 2 {dx 2 + dy 2 ) (4.17) 



10 Note that the geometry includes back-reaction from the scalar condensate, therefore the condensate 
influences the fermions even when A = 0. 
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r 2 + Q! 



with f(r 
symmetries in Eqs. 



- l+ ® , <f>{r) = fj, (l — and \i = Q. Here we have used the scaling 
[|2.4p to set L = 1 and horizon r + = 1, then the temperature is 



s< 3 



Q 2 ), 



and the dimensionless temperature is 

T 



1 (3-Q 2 
4vr Q 



(4.18) 



(4.19) 



Since the scalar ip = 0, the black hole charge Q is the only tunable parameter. 

The discussion in [9] concentrated on the spectral density of the fermions at zero 
temperature, finding a signal of a Fermi surface by the presence of a delta function 
peak in the spectral function (or equivalently, a pole in the Green's function), at critical 
value of the momentum {k = kp) and zero frequency. As expected, we find that at 
finite temperature, the pole in the Green function moves off the real axis, to complex 
values of u. In the following, we investigate how the pole evolves in the complex u 
plane as the momentum k changes. 



G22 




Figure 2: The left figure shows the G22 pole location in the complex frequency plane as func- 
tion of momentum, at different temperatures, T e ff w (0.00016,0.0016,0.008) respectively 
from top to bottom. The right figure shows a plot of ImG^ vs. Re a; at T e ff w 0.00016 for a 
range of values of k: from 0.8 to 1 for the curves from yellow to red. The curve with highest 
Im G22 has k = 0.9. 



In Fig [2] we plot the typical trajectory of the pole of the spectral function (Iu^G^)) 
on the complex oj plane for a range of temperatures, in the left panel. The position 
of the pole changes as the momentum k is varied. The Fermi momentum corresponds 
to the point of closest approach to the real axis. The right panel shows the frequency 
dependence of the spectral function for one choice of the temperature. 

From this data we can determine how the minimum distance of the pole T from the 
real axis changes with temperatur qjj. We find that T grows linearly with temperature 



(Fig 3(a) ). This is consistent with the results of and extends them beyond the low 



temperature regime. 



In all the plots of the poles in the complex frequency plane, we concentrate on the vicinity of primary 
Fermi surface, which is the range of momenta where the poles of the spectral function are closest to the real 
axis. Other branches of the curves shown exist in some cases. 
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We can also measure how the imaginary part of the pole position changes as a func- 
tion of k — kF, near the Fermi momentum. For the RN black hole at finite temperature, 
r ~ {k — kp) z , with z ~ 2. The critical exponent z does not depend strongly on the 
temperature, and is close to the value of the exponent found at [9] for T = 0. The 
dependence of V on the temperature and on k — kp is summarized in Fig |3(b)| 



G22 



G22 




Im(w-«j„) 




(a) r vs. temperature 



(b) Pole position vs. momentum 



Figure 3: Variation of the pole position with temperature and momentum. The plots are for 

G*22- 



5 Non Extremal Hairy Black Hole (T < T c ) 

We now consider the backreaction of the scalar on the black hole geometry, at tempera- 
tures below the critical one. The fermions are still treated as probes. At first we do not 
include any coupling between the scalar and fermions, concentrating on the behavior 
as a function of temperature and the fermion charge. The influence of the scalar on 
the behavior of the fermions comes through the backreaction on the geometry. We 
then turn on the quartic coupling A discussed above, and demonstrate the influence of 
changing the quartic coupling, the fermion charge and the temperature on features of 
the spectral functions of the fermions. 

Fig. U] shows a plot of Iu^G^) vs. frequency ui for zero coupling. We see the typical 
"peak-dip-hump" found in |16| . 

In Fig. [5] we plot the location of the pole in the spectral density \mG22 as function 
of momentum and frequency, for various values of the coupling, including A = 0. Note 
that unlike the high temperature case, discussed in the previous section, the spectral 
density ImGn exhibits some non-analytic behavior, though for any momentum and 
frequency the peak of ImGn is much broader than the one of ImG^- For that reason 
we concentrate exclusively on the spectral density \vnG22 below. 

The qualitative features of Fig. [5] correspond to the existence of a fairly sharp 
Fermi surface (smoothed out by finite temperature), with gapless fermionic excita- 
tions. Unlike the high temperature case, and similar to the discussion in [31 J (for zero 
temperature), we find a line of sharp excitations in the frequency- momentum plane, 
rather than an isolated point. Note that we find sharp excitations for both signs of the 
momentum, or in other words both for particle and hole excitations. 




20 40 60 80 100 



Figure 4: Im(Cr22) vs. frequency u at the Fermi momentum. This is the so-called "peak-dip- 
hump" behaviour, with some ripples in the large u region. The data is for T e ff = 0.004, qp = 
1,A = 0. 




Figure 5: Peak position {k, u) as a function of coupling. A = 0, 0.5, 1, 2 from bottom to top 
(green to red). 
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Figure 6: Pole position (k, oj) as a function of temperature. The temperature T e ff = 
0.023,0.014,0.004,0.002,0.0013,0.0009 from yellow to red (left to right). The other param- 
eters are qp — 1, A = 0. 
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Fig. [6] is similar to Fig. [5j except that here the temperature is varied for a fixed 
coupling. We can also determine the behavior of the imaginary part of the pole position 
in the complex oj plane as a function of temperature, similar to what we did in Sec. 
HI We find that in this case the temperature dependence is no longer linear, rather we 
find that for A = 0, Im u p ~ T v , with r/ ~ 5. 

Turning on and varying the scalar-fermion coupling (in addition to the fermion 
charge) corresponds to scanning different boundary theories, as opposed to variation 
of physical parameters such as the temperature. In all those boundary theories, we 
find that the qualitative properties of the system do not change much: in the con- 
densed phase we still have stable fermionic excitations (as zero temperature) which are 
gaplesO It is surprising to find such stable excitations for a large set of theories, par- 
ticularly for both signs of the coupling A. Perhaps this is somehow related to another 
feature that sets these theories apart, namely having a dual description in terms of 
classical gravity. 






Re(m) 



Figure 7: Contour plots for varying A for temperature T e ff ~ 0.004. 

In Fig [7] we plot the contour of poles for changing A, at a fixed temperature. The 
coupling A here determines the range of values of k which give rise to poles, and there 
are no gapped excitations in this picture. In the zero temperature limit , the gapless 
excitations become stable. We find that turning on A (for a fixed temperature) seems 
to stabilize the fermionic excitations, for either sign of A. For example, the critical 
exponent z defined above and in |9] increases monotonically with |A|. 

We also checked the dependence of these feature on the fermion charge qF- These 
results are summarized in figure In this case as well, we find that while qualitative 
features do not change as the fermion charge is varied, the detailed features do. 



12 These observations are consistent with the expectations in [5Tj . 
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Figure 8: The pole position as function of fermion charge qp- The plots are for T e ff = 0.004. 
As qF increases, the curve in the complex uj plane moves closer to the real axis, and the 
fermion excitations become more stable. In the left panel, lm(uj) is plotted as function as 
qF, while the right panel shows the corresponding Fermi momenta &f. 
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